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Abstract 

In the mean-field limit the dynamics of a quantum Bose gas is described by a Hartree equation. 
We present a simple method for proving the convergence of the microscopic quantum dynamics to 
the Hartree dynamics when the number of particles becomes large and the strength of the two-body 
potential tends to like the inverse of the particle number. Our method is applicable for a class of 
singular interaction potentials including the Coulomb potential. We prove and state our main result 
for the Heiscnberg-picturc dynamics of "obscrvablcs" . thus avoiding the use of coherent states. Our 
formulation shows that the mean-field limit is a "semi-classical" limit. 

1 Introduction 

Whenever many particles interact by means of weak two-body potentials, one expects that the potential 
felt by any one particle is given by an average potential generated by the particle density. In this mean- 
field regime, one hopes to find that the emerging dynamics is simpler and less encumbered by tedious 
microscopic information than the original iV-body dynamics. 

The mathematical study of such problems has quite a long history. In the context of classical 
mechanics, where the mean-field limit is described by the Vlasov equation, the problem was successfully 
studied by Braun and Hepp [3], as well as Neunzert [16]. The mean- field limit of quantum Bose gases 
was first addressed in the seminal paper [10] of Hepp. We refer to [6] for a short discussion of some 
subsequent results. The case with a Coulomb interaction potential was treated by Erdos and Yau 
in [6]. Recently, Rodnianski and Schlein [21] have derived explicit estimates for the rate of converge to 
the mean-field limit, using the methods of [10] and [Sj. A sharper bound on the rate of convergence in 
the case of a sufficiently regular interaction potential was derived by Schlein and Erdos [22], by using a 
new method inspired by Lieb-Robinson inequalities. In [151 17]. the mean-field limit (N — > oo) and the 
classical limit were studied simultaneously. A conceptually quite novel approach to studying mean- 
field limits was introduced in [8]. In that paper, the time evolution of quantum and corresponding 
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"classical" observables is studied in the Heisenberg picture, and it is shown that "time evolution 
commutes with quantization" up to terms that tend to in the mean- field ("classical") limit, which 
is a Egorov-type result. 

In this paper we present a new, simpler way of handling singular interaction potentials. It yields 
a Egorov-type formulation of convergence to the mean- field limit, thus obviating the need to consider 
particular (traditionally coherent) states as initial conditions. Another, technical, advantage of our 
method is that it requires no regularity (traditionally H 1 - or i/ 2 -regularity) when applied to coherent 
states. 

Such kinds of results were first obtained by Egorov [5j for the semi-classical limit of a quantum 
system. Roughly, the statement is that time-evolution commutes with quantization in the semi- 
classical limit. We sketch this in a simple example: Let us start with a classical Hamiltonian system 
of a finite number / of degrees of freedom. The classical algebra of observables 21 is given by (some 
subalgebra of) the Abelian algebra of smooth functions on the phase space T := M 2 ? . Let H S 21 be 
a Hamilton function. Together with the symplectic structure on T, H generates a symplectic flow (j) 1 
on r. Now we define a quantization map : 21 — > 21, where 21 is some subalgebra of B {L 2 (Rf ) ) . For 
concreteness, let (-)~ be Weyl quantization with deformation parameter K. This implies that 

[A h ,B R ] = j{A\B} h + 0(h 2 ), 

for H — > 0. The quantized Hamilton function defines a 1-parameter group of automorphisms on 21 
through 

A ^ e itSh / h Ae- itSh/h , AG 21. 
A Egorov-type semi-classical result states that, for all A € 21 and t € K, 

(A7p) n = j tS */ h A n e- itS */ h + R h (t), 

where ||-R B (t)|| -> as h -> 0. 

This approach identifies the semi-classical limit as the converse of quantization. In a similar 
fashion, we identify the mean-field limit as the converse of "second quantization". In this case the 
deformation parameter is not H, but A" 1 , a parameter proportional to the coupling constant. We 
consider the mean-field dynamics (given by the Hartree equation in the case of bosons), and view it as 
the Hamiltonian dynamics of a classical Hamiltonian system. We show that its quantization describes 
A-body quantum mechanics, and that the "semi-classical" limit corresponding to A" 1 — > takes us 
back to the Hartree dynamics. 

We sketch the key ideas behind our strategy. 

(1) Use the Schwinger-Dyson expansion to construct the Heisenberg-picture dynamics of p-particle 
operators 

e [tHN A N {a^)e~ itHN 

(in the notation of Section [3|). 

(2) Use Kato smoothing plus combinatorial estimates (counting of graphs) to prove convergence of 
the Schwinger-Dyson expansion on A-particle Hilbert space, uniformly in A and for small \t\. 
Diagrams containing I loops yield a contribution of order A~'. 

(3) Use Kato smoothing plus combinatorial estimates to prove convergence of the iterative solution 
of the Hartree equation, for small \t\. 
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(4) Show that the Wick quantization of the series in (3) is equal to the series of tree diagrams in 
(2). 

(5) Extend (2) and (3) to arbitrary times by using unitarity and conservation laws. 

This paper is organized as follows. In Section [2] we show that the classical Newtonian mechanics 
of point particles is the second quantization of Vlasov theory, the latter being the mean-field (or 
"classical") limit of the former. The bulk of the paper is devoted to a rigorous analysis of the mean- 
field limit of Bose gases. In section[3]we recall some important concepts of quantum many-body theory 
and introduce a general formalism which is convenient when dealing with quantum gases. Section 0] 
contains an implementation of step (1) above. The convergence of the Schwinger-Dyson series for 
bounded interaction potentials is briefly discussed in Section [5l Section [6] implements step (2) above. 
Steps (3), (4) and (5) are implemented in Section [71 Finally, Section [8] extends our results to more 
general interaction potentials as well as nonvanishing external potentials. 

Acknowledgments. We thank W. De Roeck, S. Graffi and A. Pizzo for useful discussions and encour- 
agement. We would also like to thank a referee for pointing out Ref. [TJ] in connection with the remark 
following Corollary 17. 4[ 



2 Mean-field limit in classical mechanics 



In this section we consider the example of classical Newtonian mechanics to illustrate how the atomistic 
constitution of matter arises by quantization of a continuum theory. The aim of this section is to give 
a brief and nonrigorous overview of some ideas that we shall develop in the context of quantum Bose 
gases, in full detail, in the following sections. 

A classical gas is described as a continuous medium whose state is given by a nonnegative mass 
density d/i(x,t>) = Mf(x,v) dxdv on the "one-particle" phase space M 3 x R 3 . Here M is the mass of 
one "mole" of gas; n{A) is the mass of gas in the phase space volume A C M 3 x ]R 3 . Let f dx dv f(x, v) = 
v < oo denote the number of "moles" of the gas, so that the total mass of the gas is /i(M 3 x R 3 ) = uM. 
An example of an equation of motion for f(x,v) is the Vlasov equation 

dtf t (x,v) = -{v-V x f t )(x,v) + -(VV eS [f t ]-V v f t )(x,v), (2.1) 

where m is a constant with the dimension of a mass, t denotes time, and 

V eS [f){x) = V(x) + J dyW(x-y) j dv f(y,v) . 

Here V is the potential of external forces acting on the gas and W is a (two-body) potential describing 
self-interactions of the gas. 

The Vlasov equation arises as the mean-field limit of a classical Hamiltonian system of n point 
particles of mass m, with trajectories (xi(t))f =1 , moving in an external potential V and interacting 
through two-body forces with potential N~ 1 W(xi — Xj). Here N is the inverse coupling constant. 
We interpret N as "Avogadro's number", i.e. as the number of particles per "mole" of gas. Thus, 
M = mN and n = vN . More precisely, it is well-known (see j3j[l6]) that, under some technical 
assumptions on V and W, 

n 

f t (x, v) = w*-lim - V 5(x - Xi (t)) 5(v - Xi (t)) (2.2) 

n—*oo n *■ — ' 
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exists for all times t and is the (unique) solution of (|2.ip . provided that this holds at time t = 0. Here, 
ft is viewed as an element of the dual space of continuous bounded functions. 

Note that n and N are, a priori, unrelated objects. While n is the number of particles in the 
classical Hamiltonian system, iV -1 is by definition the coupling constant. The mean-field limit is the 
limit n — > oo while keeping n oc N; the proportionality constant is v. 

It is of interest to note that the Vlasov dynamics (j 2 . 1 j) may be interpreted as a Hamiltonian 
dynamics on an infinite-dimensional affine phase space Tviasov To see this, we write 

f(x,v) = a(x, v)a(x, v) , 

where a(x, v), a(x, v ) are complex coordinates on IViasov For our purposes it is enough to say that 
Tviasov is some dense subspace of L 2 (M 6 ) (typically a weighted Sobolev space of index 1). On Tyiasov 
we define a symplectic form through 



ui = i / dx dv da(x, v) A da(x, v) . 

This yields a Poisson bracket which reads 

{a(x,v),a(y,w)} = {a{x,v),a(y,w)} = 0, 
{a(x,v),a(y,w)} = i5(x — y)5(v — w) . 

A Hamilton function H is defined on Tyiasov through 



(2.3) 



H(a) :=ijdxdv^,v) 



-v ■ V x H W(x) ■ V 

m 



H / dx dv a(x, v) 

m 



a(x, v) 

dy dw VW(x — y) \a(y, w)\ 2 



V v a(x,v). (2.4) 



Note that H is invariant under gauge transformations a i— > e~ ie a, a i— > e ie a, which by Noether's 
theorem implies that Jlapdxdf = J fdxdv is conserved. 

Let us abbreviate K := —W/m and F := — VW/m. After a short calculation using (|2.3p we find 
that the Hamiltonian equation of motion ctt(x,v) = {H,att(x,v)} reads 



a t (x,v) = (-v ■ V x - K(x) ■ V v )a t (x,v) - j dydw F(x - y) \a t (y,w)\ 2 ■ V v a t (x,v) 

+ dydwF(x-y)a t (y,w)a t (x,v)-V w a t (y,w). (2.5) 



Also, at satisfies the complex conjugate equation. Therefore, 
d 



dt 



a t (x,v)\ 2 = (-v ■ V x - K(x) ■ V v )\a t (x,v)\ 2 - / dydw F(x - y)\a t (y,w)\ 2 -V v \a t (x,v)\ 



We assume that 



+ \a t (x,v)\ 2 / dydw F(x - y) ■ [a t (y , w)V w a t (y , w) + a t (y , w)V w a t (y , w)] . (2.6) 



(2.7) 



|a(x,f)| = o(|(x,f)| 1 ), (x,v) — > oo . 



We shall shortly see that this property is preserved under time-evolution. By integration by parts, 
we see that the second line of (|2.6p vanishes, and we recover the Vlasov equation of motion (|2.ip for 
f = H 2 . 
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We comment briefly on the existence and uniqueness of solutions to the Hamiltonian equation of 
motion (|2.5|) . Following Braun and Hepp [3] , we assume that K and F are bounded and continuously 
differentiable with bounded derivatives. We use polar coordinates 

a = (3e lLp , 

where and /3 ^ . Then the Hamiltonian equation of motion (|2.5p reads 

Pt(x,v) = (-v-V x - K(x)-V v )p t (x,v)- J dydw F(x-y)p 2 (y,w) ■ V v f3 t (x,v) (2.8a) 
ipt(x,v) = (-v ■ V x - K(x) ■ V v )(ft{x,v) - / dydw F(x - y) $(y,w) - V v (pt(x,v) 



+ j dydw F{x-y)0 2 {y,w) ■ V w ip t {y,w). (2.8b) 
We consider two cases. 

(i) ip = 0. In this case a = (3 and the equations of motion (|2.8p are equivalent to the Vlasov equation 
for / = . The results of [3l[16] then yield a global well-posedness result. 

(ii) ip 0. The equation of motion (|2.8ap is independent of (p. Case (i) implies that it has a unique 
global solution. In order to solve the linear equation (|2.8bj) . we apply a contraction mapping 
argument. Consider the space X := {(p G C(K 6 ) : V(p £ L°°(1R 6 )}. Using Sobolev inequalities 
one finds that X, equipped with the norm ||<^||x := \<p{0)\ + HV^Hoo, is a Banach space. We 
rewrite (|2TFb|) as an integral equation, and using standard methods show that, for small times, 
it has a unique solution. Using conservation of J dxdv 01 we iterate this procedure to find a 
global solution. We omit further details. 

Note that, as shown in [3], the solution (5t can be written using a flow (ff on the one-particle phase 
space: f3 t (x,v) = /?o(<^ _ *(x, v )). The flow 0*(x,u) = (x(t),v(t)) satisfies 

x(t) = v(t), 

v(t) = K(x(t))+ [ dydw (3 2 t {y,w)F{x(t) - y) . 



Using conservation of f dxdv 01 we find that there is a constant C such that \<fr~ t (x,v)\ ^ |(x,v)|(l + 
t) + C(l + t 1 ) . Therefore (|2.7p holds for all times t provided that it holds at time t = 0. 

The Hamiltonian formulation of Vlasov dynamics can serve as a starting point to recover the 
atomistic Hamiltonian mechanics of point particles by quantization: Replace 

a(x,v) — ► a* N (x,v), a(x,v) — > ajy(x,v), 

where and a at are creation and annihilation operators acting on the bosonic Fock space T+ (L' 
see Appendix |Aj They satisfy the canonical commutation relations (|A.2[) ; explicitly, 



[a N (x,v),a N (y,w)] = [a* N (x,v),a* N (y,w)] = 0, 

1 

N 



[a N (x,v),a* N {y,w)] = — S(x - y)S(v - w) . (2.9) 
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Given a function A on IViasov which is a polynomial in a and a, we define an operator An on T+ by 
replacing a* with a* and Wick-ordering the resulting expression. We denote this quantization map 
by (-)jv- Here, iV -1 is the deformation parameter of the quantization: We find that 



[A N ,B 



N 



-i 



-{A,B} N + 0(N~ 2 ), 



for N — > oo. Here A and B are polynomial functions on IViasov 

The dynamics of a state $ G J- is given by the Schrodinger equation 



t ■ 



(2.10) 



where Hjy is the quantization of the Vlasov Hamiltonian H. In order to identify the dynamics given by 
(|2.1(jp with the classical dynamics of point particles, we study wave functions & n '(xi } vi, . . . ,x n ,v n ) in 
the n-particle sector of J-+, and interpret p^ n > := |$| 2 as a probability density on the n-body classical 
phase space. If € F + denotes the vacuum vector annihilated by aj\r(x,v) then 



N n / 2 f , , 

—j=r / dxi dv\ • • • dx n dv n & n> {xi,vi,...,x n ,v n )a N (x n ,v n )---a N (xi,vi)Q. 
vn! J 



It is a simple matter to check that (|2.9|) and (|2.10p imply that 



E 

i=l 



-Vi -V Xi + — VV( Xi ) ■ V Vi 
m 



A" 1 — ' m 



(n) 

Also, satisfies the same equation. Therefore, 



E 

i=i 



-«i • V Xi + — W(x;) • V l 
m 



Pi 



+ 1 E lv ^( 

AT ^— ' m 



(n) 



This is the Liouville equation corresponding to the Hamiltonian equations of motion of n classical 
point particles, 



d t Xi 
md t Vi 



Vi 



w ^)-atE w ^ 



Analogous results can be proven if a* N and are chosen to be fermionic creation and annihilation 
operators obeying the canonical anti-commutation relations and acting on the fermionic Fock space 
^_(L 2 (M 6 )). " 



3 Quantum gases: the setup 

Although our main results are restricted to bosons, all of the following rather general formalism remains 
unchanged for fermions. We therefore consider both bosonic and fermionic statistics throughout 
Sections [3] - [H Details on systems of fermions will appear elsewhere. 
Throughout the following we consider the one-particle Hilbert space 

H:= L 2 (R 3 ,dx) . 



G 



We refer the reader to Appendix [A] for our choice of notation and a short discussion of many-body 
quantum mechanics. 

In the following a central role is played by the p-particles operators, i.e. closed operators 
on T-l± = P±TL® P ', where P+ and P_ denote symmetrization and anti-symmetrization, respectively. 
When using second-quantized notation it is convenient to use the operator kernel of a^ p \ Here is 
what this means (see [T8] for details): Let y(R d ) be the usual Schwartz space of smooth functions of 
rapid decrease, and S^'fM. ) its topological dual. The nuclear theorem states that to every operator 
A on L 2 (R d ), such that the map (/, g) t—> (/ ,Ag} is separately continuous on ^(R^) x y(R. d ), there 
belongs a tempered distribution ("kernel") A G y'(M. 2d ), such that 

(f,Ag) = A(f®g). 

In the following we identify A with A. In the suggestive physicist's notation we thus have 

(f,a {p) g) = / dxi---dx p dyi---dypJ(x 1 ,...,x p )a ( - p) (x 1 ,...,Xp;y 1 ,...,yp)g(yi,...,y p ), 



where /, g G y{M. p ). It will be easy to verify that all p-particle operators that appear in the following 

satisfy the above condition; this is for instance the case for all bounded G B(H® P ). 

I s (p) 
Next, we define second quantization A^v- It maps a closed operator on Tt± to a closed operator 

on T± according to the formula 

A N (a^) := J dxi ■ ■ -dxpdyi ■ ■ ■ dy p 

i>*N(xp) ■ ■ ■ $* N {xi) a®> (xi, ...,x p ;yx,.. . , y p ) ^jv(yi) ■ ■ -^jv(%>) • (3.1) 

Here ip^ := where tp^ is the usual creation or annihilation operator; see Appendix [Al 

In order to understand the action of Ajv(a^) on Tt± \ we write 

]\jn/2 f ^ ^ 

$ {n) = —=■ dz 1 ---dz n $ ( - n \z 1 ,...,z n )rN(zn)---r N (zi)n 



and apply Ajv(a^) to the right side. By using the (anti) commutation relations (1A.2P to pull the p 
annihilation operators ipN(yi) through the n creation operators ip^(zi), and iPn(x) 0, = 0, we get the 
"first quantized" expression 

A,( a «)| M „„ = (*0^(« W »^')^. (3 . 2 ) 
n ± [0, n<p. 

This may be viewed as an alternative definition of Aj\i(a^). 

We define & as the linear span of {A N (a^) : p G N, oW G B(ft!j!?)}- Then ^ is a *-algebra 
of closable operators on T± (see Appendix |A|). We list some of its important properties, whose 
straightforward proofs we omit. 

(i) A N (a^T = A N ((a^)*) . 

(ii) If a(P) G B{7i [ l ] ) and G B{H ( £), then 

= E (^)(^A^). r 6W), (3.3) 
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where 

a (p). r 6(9) := P ± ( a (p) ® l (<? - r) )(l (p - r) ®6 (9) )P± G B(H^ +q ~ r) ) ■ (3.4) 

(iii) The operator A(a^) leaves the n-particle subspaces invariant. 

(iv) If W G B(H^) and 6 G B(W) is invertible, then 

T{b- l )l N {a^)T(b) = k N ((b- 1 ) m a ip) b® p ) , (3.5) 
where T(6) is defined on H± by 6® n . 

(v) If a® g B(wi p) ) then 



Aat(oW)|^») < (j^J ||aW||. (3.6) 

Of course, on an appropriate dense domain, (13 .31) holds for unbounded operators and b^ too. We 
introduce the notation 

[ a <P) , &(«)] r := G (p) . r 6 (9) - 6 (?) » r a (p) . (3.7) 
Note that [a ( - p \b^] = 0. Thus, 

We now move on to discuss dynamics. Take a one-particle Hamiltonian h^' = h of the form 
h = —A + u, where A is the Laplacian over R 3 and v is some real function. We denote by V the 
multiplication operator v(x). Two-body interactions are described by a real, even function w on 
M 3 . This induces a two-particle operator = W on Ti® 2 , defined as the multiplication operator 

w(x\ — X2). We define the Hamiltonian 

H N := A N (h) + ^A N (W). (3.9) 



n 



p 



Under suitable assumptions on v and w that we make precise in the following sections, one shows that 
Hn is a well-defined self-adjoint operator 
we have the "first quantized" expression 



(n) 

Hn is a well-defined self-adjoint operator on T±. It is convenient to introduce Hjy := NHjy. On 7i± 



H N \ n ( n) = I> + ^ E W * =■ R o + ^W, (3.10) 



N ^— ' iV 

i=l l^l<j^n 



in self-explanatory notation. 



4 Schwinger-Dyson expansion and loop counting 

Without loss of generality, we assume throughout the following that t ^ 0. 
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Let aW G B(W^) and w be bounded, i.e. «; G L°° 



Using the fundamental theorem of calculus 



and the fact that the unitary group (e ltH °)t&m. is strongly differentiable one finds 
e itHN A N (a {p) )e- itHN <S> {n) 

= qMH n e -isHo e itH A N ( a (p)\ g-W-ffOgis-ffog-is-ffjv <j>(n) I 



Jo 2 



A^(^ s ),A^(aJ p) ) 



jsH -isH N ^(n) 



where (-)j := r(e lt ' l )(-)r(e lth ) denotes free time evolution. As an equation between operators defined 
on J 7 ^., this reads 



UH N 



«A N (a ( P ) )e- itHir = A N (af>)+ / ds e' 

Jo 



isH N n -isH iN 
2 



A N (W s ),A N (a < f ) ) 



e isH 0e -isH N 



(4.1) 



Iteration of (|4.1|) yields the formal power series 



dt 



It is easy to see that, on H±\ the fc-th term of 114. 2 j) is bounded in norm by 



2 k 



A N (W tk ),... A N (W tl ),A N (a\ 



(Ph 



(4.2) 



(tn 2 ||^|| oc /iV) fe fn 



k\ 



N 



(4.3) 



Therefore, on TL± \ the series (|4.2|) converges in norm for all times. Furthermore, (|4.3|) implies that 
the rest term arising from the iteration of (|4.ip vanishes for k — > oo, so that (|4.2p is equal to (|4.ip . 

The mean-field limit is the limit n = i/iV — ► oo, where v > is some constant. The above estimate 
is clearly inadequate to prove statements about the mean-field limit. In order to obtain estimates 
uniform in N, more care is needed. 

To see why the above estimate is so crude, consider the commutator 



iN 



A N (W s ),A N (a ( f ) ) 



» 



l(«-P)]p ± . 



We see that most terms of the commutator vanish (namely, whenever p < i < j). Thus, for large n, 
the above estimates are highly wasteful. This can be remedied by more careful bookkeeping. We split 
the commutator into two terms: the tree terms, defined by 1 $C i ^ p and p + 1 ^ j $J n, and the loop 
terms, defined by 1 ^ % < j ^ p. All other terms vanish. This splitting can also be inferred from (|3.8p . 

The naming originates from a diagrammatic representation (see Figure liTTj) . A p-particle operator 
is represented as a wiggly vertical line to which are attached p horizontal branches on the left and 
p horizontal branches on the right. Each branch on the left represents a creation operator i/j^(xi), 
and each branch on the right an annihilation operator ipN(Ui)- The product Ajv(a^)Ajv(&^) of 
two operators is given by the sum over all possible pairings of the annihilation operators in Ajv(o^) 
with the creation operators in Ajv(&^)- Such a contraction is graphically represented as a horizontal 
line joining the corresponding branches. We consider diagrams that arise in this manner from the 
multiplication of a finite number of operators of the form Ajv(o^)- 
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» 



ip+l,s 



W, 



Figure 4.1: Two terms of the product A]\r(a[ p ^) An(W s ), represented as labelled diagrams. A tree term 
(left) produces a tree diagram. A loop term (right) produces a diagram with one loop. 

We now generalize this idea to a systematic scheme for the multiple commutators appearing in the 
Schwinger-Dyson expansion. To this end, we decompose the multiple commutator 



(iivy 



A N (W tk ),... A N (W tl ),A N (a 



into a sum of 2 k terms obtained by writing out each commutator. Each resulting term is a product 
of k + 1 second-quantized operators, which we furthermore decompose into a sum over all possible 
contractions for which r > in (13.3H (at least one contraction for each multiplication). The restriction 
r > follows from [a^ p \ b^]o = 0. This is equivalent to saying that all diagrams are connected. 

We call the resulting terms elementary. The idea is to classify all elementary terms according to 
their number of loops I. Write 



(LZV)* 
2 k 



A N (W t , 



A N (W tl ),A N (4 p) ) 



1 



£i^(<L*> w )). 



1=0 



N l 



(4.4) 



where G 



(k,i) 

t,t!,...,t k 



(o^') is a (p + k — Z)-particle operator, equal to the sum of all elementary terms with 



/ loops. It is defined through the recursion relation (on TC± +k l ') 



(ph 



\{p + k-l-l) 
/r p+k — I 

+ i' 



p+k-l—l 

E 

i=l 



W tk ,Gf-lt^lM P) ) 



[w^+k-^G^^iaM) 



<S> 1 



P± 



E 

l^i<j^p+k-l 



(4.5) 



as well as Gj 0,0) (aW) := a!f\ If I < 0, I > k, or p + k - I > n then G$2...,t h ( a(p) ) = °- The 
interpretation of the recursion relation is simple: a (k, ^)-term arises from either a (k — l,Z)-term 
without adding a loop or from a (k — 1, 1 — l)-term to which a loop is added. It is not hard to see, 
using induction on k and the definition (|4.5p . that (14. 4h holds. It is often convenient to have an explicit 
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formula for the decomposition into elementary terms: 



a=t 



where G^} h (a^) is an elementary term, and c(p, k, I) is the number of elementary terms in 




Figure 4.2: The labelled diagram corresponding to a one-loop elementary term in the commutator of 
order 4. 

In order to establish a one-to-one correspondence between elementary terms and diagrams, we 
introduce a labelling scheme for diagrams. Consider an elementary term arising from a choice of 
contractions in the multiple commutator of order k, along with its diagram. We label all vertical 
lines v with an index i v G N as follows. The vertical line of is labelled by 0. The vertical line 
of the first (i.e. innermost in the multiple commutator) interaction operator is labelled by 1, of the 
second by 2, and so on (see Figure FO]) . Conversely, every elementary term is uniquely determined 
by its labelled diagram. We consequently use a = 1, . . . , c(p, k, I) to index either elementary terms or 
labelled diagrams. 

Use the shorthand t = (ti, . . . , tfc) and define 

Gf'V p) ) ■= f dtGf t 'V P) )- (4-6) 

JA k {t) 

In summary, we have an expansion in terms of the number of loops /: 

oo k 1 

e^A w(a W)e-^ = ££^A w (Gf<V p) )), (4.7) 

k=0 l=0 

which converges in norm on n £ N, for all times t. 

5 Convergence for bounded interaction 

For a bounded interaction potential, ||w||oo < 00 j h is now straightforward to control the mean- field 
limit. 
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Lemma 5.1. We have the bound 



G\ k /\a®) < c(p,k, 



Furthermore, 



c(p,k,l) < 2 A 



(p + k — I) 1 (p + k — 1) 



(5.1) 



(5.2) 



Proof. Assume first that / = 0. Then the number of labelled diagrams is clearly given by 2 k p ■ ■ ■ (p + 
k — 1). Now if there are / loops, we may choose to add them at any / of the k steps when computing 
the multiple commutator. Furthermore, each addition of a loop produces at most p + k — I times more 
elementary terms than the addition of a tree branch. Combining these observations, we arrive at the 
claimed bound for c(p, k, I). 

Alternatively, it is a simple exercise to show the claim, with c(p, k, I) replaced by the bound (|5.2p . 
by induction on fc. □ 

Lemma 5.2. Let u > and t < (S^H^Hoo) -1 . Then, on TC± N \ the Schwinger-Dyson series (|4.7|) 
converges in norm, uniformly in N . 

Proof. Recall that p + k — I ^ n for nonvanishing Ajv(G^' (o^)) \jAn)- Using the symbol I{A} > defined 
as 1 if A is true and if A is false, we find 

oo k 



EE N 



k=0 1=0 



dt 



•(*) 



A N (G^/\a^))\ ( , N) 



< EE^ ± ^W^ } i(2|k||oot) fe 



k=0 1=0 

oo 



p+k—1 



klv p+k ~ l \\a {p) \ 



< J](8i/||u;||oot) fc (2^) p ||a 



(P)\ 



k=0 



1 



{2vY\\a^\\, 



1 — 8f Humect 

where we used that Ya=o (?) = 2 * \ and in particular (*) ^ 2 k . 



□ 



In the spirit of semi-classical expansions, we can rewrite the Schwinger-Dyson series to get a 
"1/iV-expansion" , whereby all /-loop terms add up to an operator of order 0(N~ l ). 



Lemma 5.3. Let t < (8 

^ll^lloo) ^ (iTid L G M. ThcTi we have on 7~C± 



j tHN A N (a^) 



-itH N 



1=0 k=l V 7 



where the sum converges uniformly in N. 



Proof. Instead of the full Schwinger-Dyson expansion (|4.2p . we can stop the expansion whenever L 
loops have been generated. More precisely, we iterate (|4.ip and use (13. 8p at each iteration to split the 
commutator into tree (r = 1) and loop (r = 2) terms. Whenever a term obtained in this fashion has 
accumulated L loops, we stop expanding and put it into a remainder term. Thus all fully expanded 
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terms are precisely those arising from diagrams containing up to L — 1 loops, and it is not hard to 
show that the remainder term is of order N~ L . 

In view of later applications, we also give a proof using the fully expanded Schwinger-Dyson series. 
From Lemma 15.21 we know that the sum converges on in norm, uniformly in N, and can be 

reordered as 

oo oo „ 

e**" £*(a«) e- 1 **" = E^E / d * M<$f (« (P) )) . 

1=0 k=l ■' A *(*) 

as an identity on H± N ^ . Proceeding as above we find 

oo 1 oo „ 

E^E/ ii(() *|s w ( G £"(^))| Hr ,| 

/ — L k — I * * 

<wtt w-.** ppi- w ( ? + r o « ii« w ii 

l=L k=l 



1 oo oo 

< 7WV> EE(P+ fc " Q L (8i>Hloo*)* (2^) p ||aW | 

^ ' l=L k=l 
_. oo oo 

= pv)T E X> + ^(HHU*)** (2^) P ||a« || 



i=L fc=0 



< ^^V( 8i/ IIHIoo*)S f L , ! rTXT (2z/) p ||a^| 

i vN ) L ^ L (l-8i/||io|| 00 *) i + lV ; " ' 

1 ^LlfSHI-ll^^n^^ 



(WV) L (1 - 8^||t(;|| 00 t) i + 2 

6 Convergence for Coulomb interaction 



where we used that J2kLo(P + k ) L %k ^ (i-J)^ ■ n 



In this section we consider an interaction potential of the form 

1 

w(x) = K- — - , (6.1) 
\x\ 

where k6 1. We take the one-body Hamiltonian to be 

h = -A, 

the nonrelativistic kinetic energy without external potentials. We assume this form of h and w through- 
out Sections [6] and [3 In Section [HJ we discuss some generalizations. 



6.1 Kato smoothing 

The non-relativistic dispersive nature of the free time evolution e ltA is essential for controlling singular 
potentials. It embodied by the following dispersive estimate, which is sometimes referred to as Kato's 
smoothing estimate, as it was first derived using Kato's theory of smooth perturbations; see [20^123] . 
Here we present a new, elementary proof, which yields the sharp constant and may be easily generalized 
to free Hamiltonians of the form (—A) 7 , where 1/2 < 7 < d/2 and d denotes the number of spatial 
dimensions. 
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Lemma 6.1. For d ^ 3 and ijj £ L (R d ) we /iai>e 

diUl^e^^ll 2 < ^IHI 2 - (6.2) 
More generally, for d ^ 2 and 7 satisfying 1/2 < 7 < d/2 we Ziaue 

di |||x|-^e-"(- A ) 7 V|| 2 < c d , 7 ||^|| 2 , (6.3) 



/or some constant Cd^ > 0. 

Remark 6.2. The constant in (|6.2p is sharp. Indeed, (|6.2p is saturated if ip is Gaussian, as can be seen 
by explicit calculation. 

Remark 6.3. At the endpoint 7 = 1/2 the dispersion law of the time evolution is uj{k) = \k\. Thus 
all spatial frequency components have the same propagation speed, i.e. there is no dispersion and the 
smoothing effect of the time evolution (which relies on the fast propagation of high spatial frequencies) 
vanishes. It is therefore not surprising that the endpoint 7 = 1/2 is excluded in (|6.3p . Similarly, the 
claim is false at the other endpoint 7 = d/2. This can be seen by noting that, for instance if ip is 
Gaussian, e~ lt ^~^ d/2 ip is nonzero in a neighbourhood of for small times. Since \x\~ d is not locally 
integrable, it follows that the left-hand side of (16. 3p is 00. 



Remark 6.4. It is easy to see that our proof of (|6.3I) remains valid if the power law potential v(x) 
is replaced with a potential v satisfying 



\v 2 (k)\ 



< 



where ~ denotes Fourier transformation. 

Proof of Lemma \6.1i The left-hand side of (|6.2p defines a quadratic form in ijj. By density, if we prove 
(|6.2p for all ij) 6 S, it follows that (|6.2p holds for all ip S L 2 . Let us therefore assume that tp S S. By 
monotone convergence, we have 



/ 



dt |||xr 1 e itA i/'|| 2 = lim/fa), 



where 

f(rj) := J dt \\lx\- 1 e itA iif e-l t2 = j dt , e" itA \x\~ 2 e iiA e"^ 

In order to write the scalar product in Fourier space, we recall (see e.g. |13| ) that, for < a < d, we 
have 

\x\-«(k) = 2 d / 2 ~ a [ 2 ' 



In particular, 

(2vr) d / 2 1 



r(f ) \k\ d ~ a ' 

- (Ott\ 

x\~ 2 (k) 



(d-2)\S d ~ 1 \ \k\ d ~ 2 ' 

where S d_1 denotes the unit sphere in M. d and | S rf 1 1 its surface measure. Thus (writing ifj instead of 
ip) we get 
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Using Fubini's theorem we get 

f(ri) = {d _2 )]sd ^ J # 1 dp 2 ?fa) bi _ 1 p2|d _ 2 ^2) / dte-^e^-i) 

1 1 _J_f D 2 — n 2 ) 2 

&P1&P2 ip(pi) i r- n rTp(p 2 )27r e 2 ^ yPl P2 > 

\Pl-P2\ V 27rr ? 

< n ^'iorf-ii / d ^ d ^ l^(Pi)l l^(P2)| i 1 — r z^=e-^ {p *- p2)2 

\Pi -P2\ v2vr?7 

where in the last step we used the inequality 2ab ^ a 2 + 6 2 and symmetry. This implies 

/fa) < n L_„ IHI 2 sup [ ^ , * * e -*«-*» a . 



(d- 


2)|S d - 1 | 




1 


(d- 


2)|S^ 1 | 




2tt 


(d- 


2)|S rf " 1 | 




2vr 


(d- 


2)|S d - 1 | 



(d-2)|s^i lirn ^ \ Pl -p 2 \ d ~ 2 v^v 

Let us write p 2 = Xp and A; := pi/X, where A > and p G S d_1 . Thus we get 

* m 9 2 ,Vii "^" 2 su p / dfc n-^^ e "^ (fc2 " 1)2 

(d-2)|S d A | \>o, P es d -iJ \k-p\ d VSttt? 
We do the integral over A; using polar coordinates: 

i-d-2 

k = y/ve, dk = ^—^ — dude, v G (0, oo) , e G S rf_1 , 

where de denotes the usual surface measure on S d_1 . This gives 

A 2 



av g(v) e 2j > 



where 



Z^ -2 /■ 1 f 



9(v) ■= — - — / de --= = - / de , . (6.4) 

-i \^ve — p\ a z 2 Jgd-i \e — p/^Jv\ a z 



d- 



Next, recall Newton's theorem for spherically symmetric mass distributions (see e.g. [13]): If fi is a 
spherically symmetric, finite, complex measure on M. d , then 

/ My) | x _ 1 y |d- 2 = j^2 J My) h\vK\*\} + / My) |^p=2i{iwi>w}- 

This yields 

1 f |S d - 1 |^ fl! ~ 2 ifv^l 

2 S d 1 if v > 1 . 



Thus, g is continuous and takes on its maximum value at 1. Since 



A 2 _Aif,_iy» 



_e 2 "' 

is an approximate delta-function centred at 1 it follows that 

f 00 A 2 A * ^ 2 f 00 A 2 * 4 ^2 

sup / df g(v) e ^ = lim / dv g(v) e ^ = <?(1) . 



A Jo 
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Thus, 



2vr o vr 



2 (d-2)\S d - 1 \ 1 1 " r " d-2 1 

This completes the proof of (|6.2p . 

The proof of (|6.3|) follows the proof of (|6.2[) up to (|6.4|) • The claim then follows from 

sup / de | /Vid-2-v < 00 ' 

for p G S^" 1 and 2j > 1. □ 

In order to avoid tedious discussions of operator domains in equations such as (|4.ip . we introduce 
a cutoff to make the interaction potential bounded. For e ^ set 

:= w(x)I { \ w{x) \^ e -i } , 

so that ||w £ ||oo ^ e -1 - Now (I6.2|) implies, for d = 3 and e ^ 0, 

|u; £ e iiA V|| 2 dt / ||W< A ^|| 2 dt ^ ttk 2 \\if}\\ 2 . (6.5) 

An immediate consequence is the following lemma. 
Lemma 6.5. Let $M G ^ n) . T/ien 

/ ||P^-e- WHo $Wf dt < ^!||^»W|| 2 . (6.6) 
Jm. 2 

Proof. By symmetry we may assume that = (1,2). Choose centre of mass coordinates X := 

(x\ + X2)/2 and £ = X2 — x\, set £, X3, • • • ,%n) := & n \xi, . . . ,x n ), and write 

( \\W[ 2 e- itHa ^\\ 2 dt = [ \\w e ^)e 2itA ^^\\ 2 dt, 
Jr Jr 

since Hq = — Ai — A2 = — Ax/2 — 2Ag and [Ax, w e (0] = 0- Therefore, by (|6.5p and Fubini's theorem, 
we find 



I \\Wf 2 e- itHo & n) \\ 2 dt = [ dXdx 3 ---dx n [ dtdC\w e ( y Oe 2itA ^ {n) (X^,x 3 ,..., 
Jr J J 

^ ^- j dXdXf-dXn j dt\& n XX,t,X Z ,...,X r U 



x n )\ 2 



2 

2 



TTK 



|$M|| 2 . Q 



2 

By Cauchy-Schwarz we then find that 

J* \\WZ B *W\\ ds ^ t 1 / 2 ^' HWg-e-^^WfdaV 2 < (j^Y V (n) ll- (6.7) 

By iteration, this implies that, for all elementary terms a, 

/t pt /^„2+\ fe/2 

dtL-.y dt fc ||Gg' 0(Q) ' £ (a^)^+ fc -')|| ^ f^J ||aW||||^+ fe -')||, 



(6.8) 
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where the superscript e reminds us that G t I '^'^(affl) is computed with the regularized potential w e . 
Thus one finds 

/ 2 \ k/2 

||G?^(a«)|| < c(p, M )(^) 

for all e ^ 0. 

Unfortunately, the above procedure does not recover the factor \jk\ arising from the time-integration 
over the A;-simplex A k (t), which is essential for our convergence estimates. First iterating (|6.6p and 
then using Cauchy-Schwarz yields a factor 1/^/kl, which is still not good enough. 

A solution to this problem must circumvent the highly wasteful procedure of replacing the integral 
over A k (t) with an integral over [0, t] k . The key observation is that, in the sum over all labelled 
diagrams, each diagram appears of the order of k\ times with different labellings. 



6.2 Graph counting 

In order to make the above idea precise, we make use of graphs (related to the above diagrams) to 
index terms in our expansion of the multiple commutator 



(my 

2 k 



A N (W tk ),... A N (W h ),A N (a[ p) ) .... (6.9) 



The idea is to assign to each second quantized operator a vertex v = 0, . . . , k, and to represent each 
creation and annihilation with an incident edge. A pairing of an annihilation operator with a creation 
operator is represented by joining the corresponding edges. The vertex has 2p edges and the vertices 
1, . . . , k have 4 edges. We call the vertex the root. 

The edges incident to each vertex v are labelled using a pair A = (d,i), where d = a, c is the 
direction (a stands for "annihilation" and c for "creation") and i labels edges of the same direction; 
i = 1, . . . ,p if v = and i = 1, 2 if v = 1, . . . , k. Thus, a labelled edge is of the form {(i>i, Ai), (t>2, A2)}. 
Graphs Q with such labelled edges are graphs over the vertex set V(Q) = {(v , A)}. We denote the set 
of edges of a graph Q (a set of unordered pairs of vertices in V(Q)) by E(Q). The degree of each (v, A) 
is either or 1; we call (v, A) an empty edge of v if its degree is 0. We often speak of connecting two 
empty edges, as well as removing a nonempty edge; the definitions are self-explanatory. 

We may drop the edge labelling of Q to obtain a (multi)graph Q over the vertex set {0, . . . ,k}: 
Each edge {(t>i,Ai), (i>2 5 A2)} € E(Q) gives rise to the edge {^1,^2} £ E(Q). We understand a path 
in Q to be a sequence of edges in E(Q) such that two consecutive edges are adjacent in the graph Q. 
This leads to the notions of connectedness of Q and loops in Q. 

The admissible graphs - i.e. graphs indexing a choice of pairings in the multiple commutator (|6,9p 
- are generated by the following "growth process". We start with the empty graph Qq, i.e. E{Qq) = 0. 
In a first step, we choose one or two empty edges of 1 of the same direction and connect each of them 
to an empty edge of of opposite direction. Next, we choose one or two empty edges of 2 of the same 
direction and connect each of them to an empty edge of or 1 of opposite direction. We continue in 
this manner for all vertices 3, . . . , k. We summarize some key properties of admissible graphs Q. 

(a) Q is connected. 

(b) The degree of each (v, A) is either or 1 

(c) The labelled edge {(v±,Xi), (^2,-^2)} G E(Q) only if Ai and A2 have opposite directions. 
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Figure 6.1: An admissible graph of type (p = 4, k = 7, 1 = 3). 



Property (c) implies that each graph Q has a canonical directed representative, where each edge is 
ordered from the a-label to the c-label. See Figure [6TT1 for an example of such a graph. 

We call a graph Q of type (p, k, /) whenever it is admissible and it contains / loops. We denote by 
@(p, k, I) the set of graphs of type (p, k, I). 

By definition of admissible graphs, each contraction in (|6.9p corresponds to a unique admissible 
graph. A contraction consists of at least k and at most 2k pairings. A contraction giving rise to a 
graph of type (p, k, I) has k + / pairings. The summand in (|6.9p corresponding to any given /-loop 
contraction is given by an elementary term of the form 



(6.10) 



2k ]\fk+V 

where the (p + k — /)-particle operator is of the form 

b (p+k-i) = p± . . . Wir ^ w (a b) l(fc -0) Wlr+ljr+utvr+i ■ ■ ■ W ikjk , H P± , (6.11) 

for some r = 0, . . . ,k. Indeed, the (anti)commutation relations (|A.2|) imply that each pairing produces 
a factor of 1/N. Furthermore, the creation and annihilation operators of each summand corresponding 
to any given contraction are (by definition) Wick ordered, and one readily sees that the associated 
integral kernel corresponds to an operator of the form (|6.1ip . Thus we recover the splitting (|4.4p . 
whereby G[ k ^ tk (a^) is a sum, indexed by all /-loop graphs, of elementary terms of the form (|6.1ip . 

As remarked above, we need to exploit the fact that many graphs have the same topological struc- 
ture, i.e. can be identified after some permutation of the labels {1, . . . , k} of the vertices corresponding 
to interaction operators. We therefore define an equivalence relation on the set of graphs: Q ~ Q' if 
and only if there exists a permutation a £ Sj, such that Q' = R a {Q). Here R a (G) is the graph defined 
by 

{(v 1 ,X 1 ),(v 2 ,X 2 )} eE(R a (g)) ^ {M«i),Ai),(<t(«2),A 2 )} £E(G), 

where cr(0) = 0. We call equivalence classes [Q] graph structures, and denote the set of graph structures 
of admissible graphs of type (p, k, I) by J2(p, k, I). 

Note that, in general, R a (G) need not be admissible if Q is admissible. It is convenient to increase 
i,k,l) to include all R a {G) where a € Sk and Q is admissible. In order to keep track of the 
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admissible graphs in this larger set, we introduce the symbol ig which is by definition 1 if Q £ &(p, k, I) 
is admissible and otherwise. Because R a {Q) ^ Q if o ^ id, 

\W(p,k,l)\ = k\ \£(p,k,l)\ . (6.12) 

Our goal is to find an upper bound on the number of graph structures of type (p,k,l), which is 
sharp enough to show convergence of the Schwinger-Dyson series (|4,2p . Let us start with tree graphs: 
/ = 0. In this case the number of graph structures is equal to 2 k times the number of ordered trees3 
with k + 1 vertices, whose root has at most 2p children and whose other vertices have at most 3 
children. The factor 2 k arises from the fact that each vertex v = 1, . . . , k can use either of the two 
empty edges of compatible direction to connect to its parent. We thus need some basic facts about 
ordered trees, which are covered in the following (more or less standard) combinatorial digression. 

For x, t E R and n £ N define 

x fx + nt\ 

^ = 7+«{ n ) (<U3 > 

as well as Ao(x,t) := 1. After some juggling with binomial coefficients one finds 

n 

Y,Mx,t)A n _ k (y,t) = A n {x + y,t); (6.14) 

k=0 

see |12j for details. Therefore 

^ A ni (xi,t) ■ ■ ■ An r (x r ,t) = A n (xi H \-x r ,t). (6.15) 

rei H Yn r =n 

Set 

„ im , , 1/1 + nm\ 1 / nm\ 

C™ := A n (l,m) = — = — — , 6.16 

1 + nvn \ n J n{m — 1) + 1 \ n / 

the n'th m-ary Catalan number. Thus we have 

V c m c m = r ( r + nm \ (617) 

ni nr r + nm\ n J y J 

\-n r =n x ' 

In particular, 

E C Z---CZ = C™. (6.18) 

niH Yn m =n-\ 

Define an m-tree to be an ordered tree such that each vertex has at most m children. The number of 
m-trees with n vertices is equal to C™ . This follows immediately from Cq 1 = 1 and from (|6.18p , which 
expresses that all trees of order n are obtained by adding m (possibly empty) subtrees of combined 
order n — 1 to the root. 

We may now compute |=2(p, fe, 0)|. Since the root of the tree has at most 2p children, we may 
express |=S(p, k, 0)| as the number of ordered forests comprising 2p (possibly empty) 3-trees whose 
combined order is equal to k. Therefore, by (|6.17p . 

\3(p,k,0)\ =2* £ Cl...C* 2p = 2 fc ^( 2p ; 3 "). (6.19) 

L„„ ^fc P \ / 



lip 



Next, we extend this result to all values of / in the form of an upper bound on |J2(p, k, l)\. 



An ordered tree is a rooted tree in which the children of each vertex are ordered. 
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Lemma 6.6. Let p,k,l G N. Then 



\J2(p,k,l)\ < 2*(*) (^ P+ k 3h ) (p + k-l) 1 . (6.20) 

Proof. The idea is to remove edges from Q G fc, I) to obtain a tree graph, and then use the special 
case (pU9j) . 

In addition to the properties (a) - (c) above, we need the following property of &(p, k, I): 

(d) If Q G &(p,k,l) then there exists a subset V C {1, . . . , A;} of size I and a choice of direction 
5 : V — > {a, c} such that, for each u G V, both edges of w with direction are nonempty. 
Denote by £{v) C -E(^) the set consisting of the two above edges. We additionally require 
that removing one of the two edges of £{v) from Q, for each v G V, yields a tree graph, with 
the property that, for each v G V, the remaining edge of £{v) is contained in the unique path 
connecting v to the root. 

This is an immediate consequence of the growth process for admissible graphs. The set V corresponds 
to the set of vertices whose addition produces two edges. Note that property (d) is independent of the 
representative and consequently holds also for non-admissible Q G &(p,k,l). 

Before coming to our main argument, we note that a tree graph T G @(p, k, 0) gives rise to a 
natural lexicographical order on the vertex set {1, . . . , k}. Let v G {1, . . . , k}. There is a unique path 
that connects v to the root. Denote by = vi, V2, ■ ■ ■ , v q = v the sequence of vertices along this path. 
For each j = 1, . . . ,q — 1, let Xj be the label of the edge {vj,Vj+i} at vj. We assign to v the string 
S(v) := (Ai, . . . , Aq_i). Choose some (fixed) ordering of the sets of labels {A}, for each v. Then the 
set of vertices {1, . . . , k} is ordered according to the lexicographical order of the string S(y). 

We now start removing loops from a given graph Q G &(p, k, I). Define 1Z\ as the graph obtained 
from Q by removing all edges in \J ve y£(v). By property (d) above, TZ\ is a forest comprising I trees. 
Define T\ as the connected component of 1Z\ containing the root. Now we claim that there is at least 
one v G V such that both edges of £(v) are incident to a vertex of T\. Indeed, were this not the case, we 
could choose for each v G V an edge in £ (v) that is not incident to any vertex of T\. Call TZ\ the graph 
obtained by adding all such edges to 1Z\. Now, since no vertex in V is in the connected component 
of 72-i, it follows that no vertex in V is in the connected component TZ^. This is a contradiction to 
property (d) which requires that TZ^ should be a (connected) tree. 

Let us therefore consider the set V of all v G V such that both edges of £(v) are incident to a 
vertex of T\. We have shown that V ^ 0. For each choice of v and e, where v G V and e G £(v), we 
get a forest of I — 1 trees by adding e to the edge set of 1Z\. Then v is in the same tree as the root, 
so that each such choice of v and e yields a string S(v) as described above. We choose v\ and e{v\) 
as the unique couple that yields the smallest string (note that different choices have different strings). 
Finally, set Q\ equal to Q from which e{v\) has been removed, and Vi := V \ {v}. 

We have thus obtained an (I — l)-loop graph Q\ and a set Vi of size I — 1, which together satisfy the 
property (d). We may therefore repeat the above procedure. In this manner we obtain the sequences 
vi,...,vi and Gi,-..,Gi- Note that Qi is obtained by removing the edges e(v±), . . . , e(vi) from Q, 
and is consequently a tree graph. Also, by construction, the sequence v\, . . . ,V[ is increasing in the 
lexicographical order of Q\. 

Next, consider the tree graph Q[. Each edge e(vj) connects the single empty edge of Vj with 
direction S(vj) with an empty edge of opposite direction of a vertex v, where v is smaller than Vj in 
the lexicographical order of Q\. It is easy to see that, for each j, there are at most (p + k — I) such 
connections. 
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We have thus shown that we can obtain any Q S &(p, k, I) by choosing some tree Q\ G @(p, k, 0), 
choosing I elements Vj out of {1, . . . , k}, ordering them lexicographically (according to the order of Q{) 
and choosing an edge out of at most (p + k — I) possibilities for v±, . . . ,vi. Thus, 



\&(p,k,l)\ < (f)(p + k-l) l \9(j>,k,0)\. 
The claim then follows from (|6.12p and (16. 19h . □ 

6.3 Proof of convergence 

We are now armed with everything we need in order to estimate /^(t) d£ (affl). Recall that 

Gtt.,M P) ) = % E ^a(« W ), (6-21) 

where (a^) is an elementary term of the form (|6,lip indexed by the graph Q. We rewrite 

this using graph structures. Pick some choice V : £?(p,k,l) — ► &(p,k,l) of representatives. Then we 
get 

Qe3(p,k,l) (reS k 

Now, by definition of R a , we see that 

(k,l)(R<r(g)) (n (p)^ = n (k,l)(g) (a (p)) 



M CT (1) , — ,t(7(k) 



Thus, 



2 26i2(p,*,!)«S* " "' 

I E / a , (t /i« ( T"^»). 



Q6.2(p,fc,0 " "cW 



where 



AqW := {(h,...,t k ) :3a£S k : i RAV{Q)) = 1, (t a{1) , . . . ,t a(k) ) £ A k (t)} C [0,*]' 
is a union of disjoint simplices. 
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Therefore, (gETJ and (loTTTj) imply, for any e ft^ +fc /} , that 



1 ^ /vr K 2 A fe/2 , 



■0l 



Q6.3(p,fc,0 

« ( 2p t 3t ) (*)&-+*-')'(^)' / Vwiin* < ' +i -'»ii, 

where the last inequality follows from Lemma 16.61 Of course, the above treatment remains valid for 
regularized potentials. We summarize: 

|| <«>, ((l(rt) || ^ (* + »)(J)(p + (6.22) 

for all e ^ 0. 

Using (|6.22p we may now proceed exactly as in the case of a bounded interaction potential. Let 

P{K ^ ):= U8^- (6 - 23) 
The removal of the cutoff and summary of the results are contained in 

Lemma 6.7. Let t < p(pz, v). Then we have on TC± 

oo k -. 

e*A w (a(ri)e-^ = ££ ^ A*(G?^(a«)) , (6.24) 

k=0 1=0 

in operator norm, uniformly in N . Furthermore, for L 6 N, we have the 1/N -expansion 

1=0 k=l ^ ' 

where the sum converges on uniformly in N. 

Proof. Using (|6.22p we may repeat the proof of Lemma l5.3l to the letter to prove the statements about 
convergence. Thus (j6.24|) holds for all e > 0. The proof of (I6.24|) for s = follows by approximation 
and is banished to Appendix [Bj □ 

7 The Mean-Field Limit 

In this section we identify the mean-field dynamics as the dynamics given by the Hartree equation. 
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7.1 The Hartree equation 

The Hartree equation reads 

idftp = hip + (w* \ip\ 2 )ip. (7.1) 

It is the equation of motion of a classical Hamiltonian system with phase space T := ff 1 (lR 3 ). Here 

i^QR 3 ) is the usual Sobolev space of index one. In analogy to Ajv we define A as the map from closed 

(p) 

operators on 7i + to functions on phase space, through 
A(a (p) )(-0) := (^ p ,a (p) ^ p ) 

dxi ■ ■ ■ dx p dyi ■ ■ ■ dy p 4>(x p ) ■ ■ ■ $(xi) a (p) (xi, . . . ,x p ;yi, . . . ,y p ) ip{y{) ■ ■ ■ ip(y p ) . 



We define the space of "observables" 21 as the linear hull of {A(aW) : p G N, oW g 
The Hamilton function is given by 

H := A(h) + l -k(W) , 

i.e. 

H(rP) = Jdx\V^\ 2 + ^Jdx(w*\^\ 2 M 2 = (^,h^) + ^ 2 ,W^ 2 ). (7.2) 

Using the Hardy-Littlewood-Sobolev and Sobolev inequalities (see e.g. [15] ) one sees that is 
well-defined on T: 



l^)| 2 |W! 2 < 111,121,2 



4 s ll„/.l|4 



where the symbol < means the left side is bounded by the right side multiplied by a positive constant 
that is independent of ift. 

The Hartree equation is equivalent to 

id t i> = dfHty) . 

The symplectic form on T is given by 

u = i J dx dip(x) A dip(x) , 

which induces a Poisson bracket given by 

bl>{x),${y)} = iS(x-y), {ip(x),tl}(y)} = {$(x), ip{y)} = 0. 
For A, B G 21 we have that 

{A,B}(ip) = i J dx [^)fl^(^)-^flW8^W] . 

The "mass" function 

JV(^) := / dx |V>| 2 



is the generator of the gauge transformations ^ i— > e~ ie ip. By the gauge invariance of the Hamiltonian, 
{H, N} = 0, we conclude, at least formally, that N is a conserved quantity. Similarly, the energy H 
is formally conserved. 

The space of observables 21 has the following properties. 
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(i) A(aW) = A((oW)*). 

(ii) If aW G } ) and 6 G then 

A(oW)(6^) = A((b*)® p a^b® p )(il>) . 

(iii) If and are p- and g-particle operators, respectively, then 

(A(o (p) ),A(6^)} = ipgAQaW.&^J . (7.3) 

(iv) If a W G then 

|A(a«)ty)| < ||a {p) ||||V'l| 2p . (7.4) 

The free time evolution 

is the Hamiltonian flow corresponding to the free Hamilton function A(h). We abbreviate the free 
time evolution of observables A G 21 by A t := Ao (fa. Thus, A(a^) t = A(a| p ' ) ). 

In order to define the Hamiltonian flow on all of L 2 (IR 3 ), we rewrite the Hartree equation (|7.ip 
with initial data ip(0) = tp as an integral equation 

${t) = e~ ith ip-i f dse-^-^iw*^)] 2 )^). (7.5) 



o 



Lemma 7.1. Let -0 G L 2 (M 3 ). Then f|7.5|) /ias a unique global solution ip(-) G C(R; L 2 (]R 3 )), w/iic/i 
depends continuously on the initial data ip. Furthermore, \\ip(t)\\ = \\tjj\\ for allt. Finally, we have a 
Schwinger-Dyson expansion for observables: Let affl G v > and t < p(k,v). Then 

oo 

A(a«)Mt)) = £ A(Gf' 0) (a (p) ))(^) 

fe=0 

oo . 

= E^fc / dt{A(W tk ),...{A(W tl ),A(af)}. ..}&), (7.6) 

, „ Z JA k (t) 



k=0 



uniformly in the ball B u := {ip G L 2 (R 3 ) : ||*0|| 2 ^ v\. 

Proof. The well-posedness of (|7.5p is a well-known result; see for instance [U[23]. The remaining 
statements follow from a "tree expansion", which also yields an existence result. We first use the 
Schwinger-Dyson expansion to construct an evolution on the space of observables. We then show that 
this evolution stems from a Hamiltonian flow that satisfies the Hartree equation (|7.5p . 

First, we generalize our class of "observables" to functions that are not gauge invariant, i.e. that 
correspond to bounded operators G B{H p + ]H q + ). We set A(a^))(y>) := , a^V® p ), and 

denote by 21 the linear hull of observables of the form A(a^ q ' p ^) with a^ q,p "> G B{7i p + ; Ti 9 + ). 

It is convenient to introduce the abbreviations 

G := {A(h),-}, D := \{A(W),-}. 
Then e Gt is well-defined on 21 through (e Gt A)(if)) = A(e~ lh i/j) : where A G 21. Note also that 

D s := e Gs De- Gs = \{A{W S ), •}. 
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Let A G 21. We use the Schwinger-Dyson series for e( G + D ) t to define the flow S(t)A through 

oo „ 

S(t)A ■= y dtD tk ---D h e Gt A 

oo 

= E / d^{A(W tfc ),...{A(W tl ),A t )}...}. 



(7.7) 



fc=0 ' 



Our first task is to show convergence of (17. 7ft for small times. 

Let ^4 = A(a( 9,p )). As with (|7.3p one finds, after short computation, that 

i{A(W),A(o(«*))} = AA^W i9+1 (a^)®l) -i^(a^) 8)1)^+1") . (7.8) 

^ i=i j=i ' 

Thus we see that the nested Poisson brackets in (|7.7p yield a "tree expansion" which may be described 



as follows. Define T f ^ t k ( a ^ q,P ^) recursively through 



T f (0) (a(^)) .- ^ , 

i=l 
p+fc-1 



a 



(<?>p) 



Note that ^§ 4(0^) is an operator from wf +/cj to H^ +k> . Moreover, (72]) implies that 



f 



2 ^{A(mj,...{A(^ tl ),A(ar ; )}---} = A(«.., t (a^) 
Also, by definition, we see that for gauge- invariant observables we have 



(7.9) 



4Uw = Ca( oW ) 



,(fe,0) 



,(p)> 



We may use the methods of Section [6] to obtain the desired estimate. One sees that t (a^) is 

a sum of elementary terms, indexed by labelled ordered trees, whose root has degree at most p + q, 
and whose other vertices have at most 3 children. From (|6.17p we find that there are 

p + q fp + q + 3k 
p + q + 3k\ k 

unlabelled trees of this kind. Proceeding exactly as in Section [6] we find that 

'p + q + 3k\ /TTK 2 t\ k/2 



dt 



where & p+k 1 G H { l +k) 



(<J,p)llll<l)(p+ fc )l 



Let ip G 1/ 



with ||t/>|| 2 ^ v. Then |A(a^)(^)| < \\a^ \\ U\\ p+q implies 



/ dt i{A(W t J,...{A(W tl ),A(a^ ) )}...}W 



< 



p + q + 3k^ ( ttkH^ k/2 ||a(w0 n ^(p+gj/a 



(7.10) 
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Convergence of the Schwinger-Dyson series (|7.7p for small times t follows immediately. 

Thus, for small times t, the flow S(t) is well-defined on 21, and it is easy to check that it satisfies 
the equation 



S{t)A = e Gt A+ [ ds S(s)De G ^- a) A, (7.11) 
Jo 



for all A G 21. 

In order to establish a link with the Hartree equation (|7.5p . we consider / G L 2 (]R 3 ) and define the 
function Fj G 21 through i ? /(V') : = (f^ip)- Clearly, the mapping / i— > (S(t)Ff)(ip) is antilinear and 
(|7.10p implies that it is bounded. Thus there exists a unique vector ip{t) such that 

{s{t)F f m = : (f,m)- 

We now proceed to show that (S(t)A)(ip) = A(ip(t)) for all A G 21. By definition, this is true for 
A = Ff. As a first step, we show that 

S(t)(AB) = (S(t)A)(S(t)B), (7.12) 

where A, B G 21. Write 

oo „ 

5(t)(AS) = ^ / dt A fc ---A 1 e G *(A5) 

oo „ 

= £ / dt A fe • • • Ai (^*#t) , 

where we used e G *(j4.B) = (e Gt A)(e Gt B). We now claim that 

I dtD tk ---D tl (A t B t ) = V / dt / ds (A t • • • A-i^t) (A m • • • D sl B t ) , (7.13) 

where the sum ranges over Z,m ^ 0. This follows easily by induction on k and using D S (AB) = 
A(D S B) + (D S A)B. Then ([715]) follows immediately. 

Next, we note that (|7.12p implies that (S(t)A)(ip) = A(ip{t)), whenever A is of the form A = 
A(a^), where 

^2 p+\fi ® • • • ® HWi ® • • • ® > ( 7 - 14 ) 



a (9,P) 



where the sum is finite, and f-,g{ G L 2 (M 3 ). Now each a^' p ^ G B(Txf ,7tf ) can be written as the 

,(?<P)l 



weak operator limit of a sequence (a^'^OneN of operators of type (|7.14p . One sees immediately that 

limA(a^ p) )(^(t)) = A(a (9 ' p) )(^(t)). 

n 

On the other hand, uniform boundedness implies that sup n ||ai 9 ' p ^ || < oo, so that 

q jMq+k) tjt . tit . fa 11 ( fc )"| W. , . . . W. . , ,/,®(p+*0\ 



C \\ai q ' p) \ 



w irJr , tvr ■ ■ ■ w iin , tvi r ig+h) w lr+ljr+utvr+i ■ ■ ■ w lk3k , tvk r {p+k) 
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justifies the use of dominated convergence in 

lim(S(t)A(a^)))(VO = (S(t)A(a(w>))(V0. 

n 

We have thus shown that 

(S(t)A)(if>) = A(i/>(t)), V-AeB. (7.15) 
Let us now return to (|7. 1 1 j) . Setting A = Ff , we find that (|7.1ip implies 

(f,m) = (f,e- ih ^)+^ds±(s(s){A(W),(F f ) t _ s })&) 
= </,e-*V)+ /*d a ({A(W0.(i ? /)*-.})^W). 

JO 

where we used (|7.15p . Using (I7.8p we thus find 

</,Y>(*)> = (/,e-^)-i / d*<(e^'-'V)®^(*),W#0®^(*)>, (7.16) 

■/ o 

which is exactly the Hartree equation (|7.5p projected onto /. We have thus shown that ift(t) as defined 
above solves the Hartree equation. 

To show norm-conservation we abbreviate F(s) ■= (w * \i/j(s)\ 2 )i/j(s) and write, using (|7.5p . 

II^)I| 2 -||^H 2 = i fas [(F(s),e- ish ^)-(e- ish ^,F(s))] 
Jo 

+ [ ds [ dr (e ish F(s),e irh F{r)) . 
Jo Jo 

The last term is equal to 

I ds [ S dr [(e ish F(s),e irh F(r)) + (e irh F{r),e ish F(s))} . 
Jo Jo 

Therefore (|7.5|) implies that 

\\i;(t)\\ 2 -U\\ 2 = i [ t ds(F(s),i;(s))-i fds^^is)) = 0, 
Jo Jo 

since (F(s) , ip(s)) 6 R, as can be seen by explicit calculation. Thus we can iterate the above existence 
result for short times to obtain a global solution. 

Furthermore, (|7.16p implies that ip(t) is weakly continuous in t. Since the norm of ip(t) is conserved, 
ip{t) is strongly continuous in t. Similarly, the Schwinger-Dyson expansion ([7.7]) implies that the map 
ip i — > ip(t) is weakly continuous for small times, uniformly in in compacts. Therefore, the map 
ip i— > ip(t) is weakly continuous for all times t, and norm-conservation implies that it is strongly 
continuous. □ 
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7.2 Wick quantization 



In order to state our main result in a general setting, we shortly discuss how the many-body quantum 
mechanics of bosons can be viewed as a deformation quantization of the (classical) Hartree theory. The 
deformation parameter (the analogue of h in the usual quantization of classical theories) is 1/N. We 
define quantization as the linear map (-) N : 21 — > 21 defined by the formal replacement ip{x) \— > iPn{x) 
and ijj(x) i— > ij)* N {x) followed by Wick ordering. In other words, 

V) N : A(aW) ~ £ N (aM). 

Extending the definition of {-) N to unbounded operators in the obvious way, we see that Hn is the 
quantization of H. 

Note that ([33]) and CL3]) imply, for A, B G 21, 

^ ^ N- 1 - 

[A N ,B N ] = —r-{A,B} N + 

so that 1/N is indeed the deformation parameter of (-)at- 




7.3 The mean-field limit: a Egorov-type result 

Let 0* denote the Hamiltonian flow of the Hartree equation on L 2 (IR 3 ). Introduce the short-hand 
notation 

a*A := Aotf, Ag21, 
S*A := e^Ae"™", AG 21. 

We may now state and prove our main result, which essentially says that, in the mean-field limit 
n = vN — > oo, time evolution and quantization commute. 

Theorem 7.2. Let A G 21, v > 0, and e > 0. T/ien i/iere exists a function A{t) G 21 stzc/i that 

sup||o!*A- A{t)\\ Loo , Bv) e, 

as we// as 

11(3*^-^)1^11 < e+ C(£ ^ M) . 

Remark. The "intermediate function" ^4(i) is required, since the full time evolution a* does not leave 
21 invariant. 

Proof. Most of the work has already been done in the previous sections. Without loss of generality 
take A = A(a^) for some p G N and a® G B(7i^). Assume that t < p(n, v). Taking L = 1 in (j6T25|) 
we get 

S 4 A/v(a^)| = £ A N (Gl k >%M))\ +o(±) . (7.17) 
Comparing this with (j7.6f) immediately yields 



3*Ajr(a«) = [a'A^)]^ + o(± 
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on TC^ N \ where [a*A(a^)] N is denned through its norm-convergent power series. This is the state- 
ment of the theorem for short times. 

The extension to all times follows from an iteration argument. We postpone the details to the 
proof of Theorem 17.31 below. In its notation A(t) is given by 

fcl=0 k m =0 

The result may also be expressed in terms of coherent states. 

Theorem 7.3. Let a® £ B(H^), ip £ L 2 (R 3 ) with \\ip\\ = 1, and T > 0. Then there exist constants 
C, (3 > 0, depending only on p, T and k, such that 

(y 9N , j tHN A N {a(P))e- itH »ip® N ^ - <^(t) 0p ,a( p ^(t) 0p >| te[0,T]. (7.18) 

Here ip(t) is the solution to the Hartree equation (|7.5p with initial data ip. 
Proof. Introduce a cutoff K 6 N and write (in self-explanatory notation) 

K-l 

N 



a T A N (a^) = ^A N (G^°\a^))+a^ K A N (a^) + ^R N>T (a^), (7.19) 



k=0 
K-l 

a T A(a^) = E A((4 fc ' 0) (a (p) )) +a^A(o (p) ). (7.20) 

k=0 

(k 0) 

To avoid cluttering the notation, from now on we drop the parentheses of the linear map GP J . We 
iterate (|7.19j) m times by applying it to its first term and get 

A'i-1 K m -1 

(STA w (a w ) = E ■■ E An^G^'^G^- 1 ^ ■ ■ ■ G^'^a^ 

ki=0 k m =0 

m-lKi-1 Kj-l 

j=l fc 1= kj=0 

+ (a T ) m ~ 1 Rn,t ( a ^ ) + ^ E E • • • E (u T ) m - l - ] RN,T (G^' 0) • • • G^-^W) . (7.21) 

j=l fci=0 fcj=0 

A similar expression without the third line holds for (a T ) m A(a^). 

In order to control this somewhat unpleasant expression, we abbreviate 



x := 



P (k,1) ' 

Assume that x < 1. Then (|6.22p and (|6.25p imply the estimates, valid on TL^\ 

|| G (_M) a (p)|| <; 4 p ||a (p) ||x fc , 



\a T >K A N (a^)\\ ^ 4 p ||a (p) " 



\RnA^ p) )\\ < (4e) p | 



1 — x 

0) II E_ 



(1-s) 



3 ' 
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Furthermore, ()7.6j) implies that 

||a^A(aW)|| Lao(fli) < na^Wj 

We also need 



,A 



jy...(jy-p + i) 

ATP 



|A(a«)M| 



N---{N-j) AT...(JV-j + l) 



ATi+1 



,(p)l 



(7.22) 



Armed with these estimates we may now complete the proof of Theorem 1 7. 3 1 Suppose that 1/2 ^ i< 1. 
Then 

A'l-l Km-1 



fci=0 fc m =0 



<C L( p + K 1 + ---+ K m f 4m(p+K 1+ -+K m ) || a (p) || _ 



Similarly, the second line of (|7.2ip on and its classical equivalent on B\ are bounded by 



x K] 4 j{p+Kl 



||„(p)| 



Finally, the last line of (|7.2ip on 7i\ is bounded by 



- V 4(i+ 1 )(p+^i+-+^i-i) || (p) I 



2V 

Now pick m large enough that T sC mr. Then it is easy to check that there exist ai, . . . , a m such 
that setting 

Kj = ctj log TV, j = 1, . . . ,m 

implies that the three above expressions are all bounded by CN'^Wa^ ||, for some f3 > 0. This remains 
of course true for all m' ^ m. Since any time t ^ T can be reached by at most m iterations with 
1/2 ^ x < 1, the claim follows. □ 

We conclude with a short discussion on density matrices. First we recall some standard results; 
see for instance [15]. Let T G £ , where £ x is the space of trace class operators on some Hilbert space. 
Equipped with the norm ||r||i := Tr|r|, C 1 is a Banach space. Its dual is equal to B, the space of 
bounded operators, and the dual pairing is given by 



(A,F) = Tr(AT), A £ B ,F G C l 
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Therefore, 



|r||i = sup |Tr(^r)|. 

Ae6, ||A||<1 



(7.23) 



Consider an iV-particle density matrix ^ Fn E C l (TL^) that satisfies TrT^ = 1 and is sym- 
metric in the sense that T^P + = T^. Define the p-particle marginals 



r 



(p) . 

N - 



TV 



p+l,...,N 1 N , 



where Tr„+i ... jy denotes the partial trace over the coordinates p + 1, . . . , N. Define furthermore 



F N (t) = e- itHN T N j tH > 
(p), 



as well as the p-particle marg inals T%>(t) of T N (t) 
Noting that 



Tr(A N (aM)T N (t)) = Tr(a«r«(t)) = Tr(a«r«(t)) + O (1) 



we see that (|7.23p and Theorem 17.31 imply the following result. 



Corollary 7.4. Let ip E 7i with ||^|| = 1, and let ip(t) be the solution of (|7.5p with initial data ip. Set 
Tn '■= (|^)(^|) ■ Then, for any p E N and T > there exist constants C, (3 > 0, depending only on 
p, T and k, such that 



Remark. Actually it is enough for IV to factorize asymptotically. If (Tn)n^ is a sequence of sym- 
metric density matrices satisfying 



lim | r 

AT— ►oo ' 



(1) 
N 



0. 



then one finds 



lim 



r$(t)-\m)(m\ 



0. 



t E 



This is a straightforward corollary of the proof of Theorem 17.31 By an argument of Lieb and Seiringer 
(see the remark after Theorem 1 in |14|). this implies that 



lim 



r^(t)-(\m)mt)\y 



0. 



t E 



for all p. 



8 Some Generalizations 

In this section we generalize our results to a larger class of interaction potentials, and allow an external 
potential. For this we need Strichartz estimates for Lorentz spaces. We start with a short summary 
of the relevant results (see [Hill] ) . 
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For 1 ^ q ^ oo and < 9 < 1 we define the real interpolation functor (•, -)g q as follows. Let Aq 
and A\ be two Banach spaces contained in some larger Banach space A. Define the real interpolation 
norm 

\[S™(t-°K{t,a)) q <it/t] ll \ q < oo, 
su Pt>0 t~ 6 K(t, a) , q = oo. 



\{A ,Ai)e, q :z 



■ 



where 



if(t,a) := inf (||ao|U + *IK |Ui) 

a=ao+ai ' 



Define (Aq, Ai)g q as the space of a £ A such that ||a||( J 4 0jJ 4. 1 ) e < oo. Then (Aq, Ai)g^ q is a Banach 
space. The Lorentz space L p,q (M? ,dx) = L p ' q is defined by interpolation as 



m := (w°,m) e , q , 



where 1 ^ Pq,P\ ^ oo, po 7^ Pi-, and 



1 1 



+ — 

Po Pi 



Lorentz spaces have the following properties that are of interest to us. First, L p ' p = L p . Second, 
jjj,oo _ jj w ^ where L v w is the weak L p space (see e.g. [Eld]). I n particular, the Coulomb potential in 
3 dimensions satisfies 

1 



G L 



3.oo 



Finally, Lorentz spaces satisfy a general Holder inequality (see p2]): Let 1 < p,p\,pi < oo and 
1 ^ q, 3i, 32 ^ oo satisfy 

111 111 
— + — = -, — + — = -. 

pi P2 p 3i 12 q 



Then we have 



We need an endpoint homogeneous Strichartz estimate proved in [TT]. For a map / : 
define the space-time norm 

1 V r 



TT JP,1 



dt \\f(t)\\ 



r 

LP>1 



Then Theorem 10.1 of [TT] implies that 



whenever 2 < r < oo and 



2 3 3 

- + - = o 
r p 2 



L 2 



We are now set for proving a generalization of (|6.2p . 

Lemma 8.1. Lei u> £ L^, + L°° . Then there is a constant C = C(w) > 0, such that 

f 

Jo 



\we 



itA M 2 dt < CU\\ 2 . 



L p ' q we 



32 



Proof. Let w = w\ + W2 with w\ G L°° and G L^. Then 



The first term is bounded by || , u;i||l°°||V'IIl 2 - To bound the second we use (18. ip and (18. 2\i with r = 2 
and p = 6 to get 



\ui2e ltA ip\\ L 2 L2 < ll^ll^ocHe 1 ^ ^|| r?r 6,2 < ||u;2||l3,oc.||'0|| j l2 . 



I L t L x 



Therefore, 



we 



AtA 



L 1 



□ 



Now let us assume that v, w G L°° + Set i?o| w («) := Y17=i Then the required generaliza 



tion of Lemma 16.51 is 

Lemma 8.2. There exists a constant C = C(w,v) such that 



f 

Jo 



|W^-e- itffo $ (n) || 2 dt < C||$ (n) || 2 



f 1 

/ \\Vie^ tHo ^\\ 2 dt < C||$( n >| 
■/ o 



w/iere $W G . 

Proof. The claim for V follows immediately from Lemma 18.11 The estimate for W follows similarly 
by using centre of mass coordinates. □ 

Finally, we briefly discuss the changes to the combinatorics arising from an external potential. We 
classify the elementary terms according to the numbers (k,l,m), where k is the order of the multiple 
commutator, I is the number of loops, and m is the number of F-operators. Thus, instead of (|4.5p . 
we have the recursive definition 



+ ip + V"lK.Cir ) (^ ) ) 



Glt U ) " ,) 1 (^)8l]P ± 



p+k— I— m— 1 

= ' lP ± E \WiP+k-l-m,t k ,^i,t x ,:.:,t k 
1=1 

li^iKji^p+k—l—m 
p+k—l—m 

+i p ± [ViA.,^f.s:^))p ± , 
i=i 

as well as "'(o^) := a^. We also set G[ k ^' m \ k (a^) = unless ^ I ^ k — m. It is again an 
easy exercise to show by induction on that 



(iN) h 
2 k 



A N (W tk ),...\A N (W tl ),A N (a ( * ) ) 



k k-l 



1=0 m =0 



EE^,(GSl(aW)). 
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Note that GjjJ^ (a^) is & p + k — I — m particle operator. 

The graphs of Section [6] have to be modified: Each vertex corresponding to a ^/-operator has one 
edge for each direction d = a,c (see Figure I87T1) . 




Figure 8.1: An admissible graph of type (p = 4, k = 7, 1 = 2, m = 2). 

Let us first consider tree graphs, 1 = 0. Take the set of trees without an external potential 
as in Section [6] By allowing each vertex v = l,...,k whose edges (a, 2) and (c, 2) are empty to 
stand for either an interaction potential W or an external potential V, we count all trees with an 
external potential. Thus, for a given m, there are at most \^(p, k, 0)| tree graphs contributing to 

G[ k ^' m \ k (a^). If / > we repeat the argument in the proof Lemma 16.6} and find that the number of 

graph structures contributing to G[ k /' m \ (a^) is bounded by 

Putting all this together, we find that 

\\G^ m \a^)\\ < QQ) ( 2p + 3A; )( P + fc-^-my(C^|| a (p)||. 

Using the condition p + fe — I — m ^ n, it is then easy to see that all convergence estimates remain 
valid with the additional factor 2 k . 

In summary, all of the results of Sections [6] and [7] hold if 

v,w G + L°° . 

A Second Quantization 

We briefly summarize the main ingredients of many-body quantum mechanics and second quantization. 
See for instance [2] for an extensive discussion. 

Let H = L 2 (R d ,dx) be the "one-particle Hilbert space", where d £ N. Many-body quantum 
mechanics is formulated on subspaces of the n-particle spaces H® n . Let = P± be the orthogonal 
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projector onto the symmetric/antisymmetric subspace of 7i.® n , i.e. 

(P ± ^)( Xl ,...,x n ) := l^(±l)H$W(x ff 

n.i *■ — » 



.(l)c--iVnV' 

crSSn 



where |<r| denotes the number of transpositions in the permutation <r, and $^ n ^ G W^ n . We define the 
bosonic n-particle space as Tt^ := P+Tl® 11 , and the fermionic n-particle space as Ti^ := P_Tt® n . 
We adopt the usual convention that = C. 
We introduce the Fock space 



oo 
n=0 



n=0 

A state <1> G JP± is a sequence <1> = ($( n ))^_ j where G T^j™ . Equipped with the scalar product 

oo 

= ^ ($("), 

n=0 

J/ 7 .)- is a Hilbert space. The vector il := (1,0,0,...) is called the vacuum. By a slight abuse of 
notation, we denote a vector of the form «3> = (0, . . . , 0, &( n \ 0, . . . ) G T± by its non- vanishing n- 
particle component $( n ). Define also the subspace of vectors with a finite particle number 

J± := {<]? G T± : <£ (n) = for all but finitely many n} . 

On T± we have the usual creation and annihilation operators, ip* and ip, which map the one- 
particle space H into densely defined closable operators on T±. For / G H and G J r ±, they are 
defined by 



n 

(^(/)cD) (n) (x 1 ,...,^) := T j](±l) i - 1 /(^ M (H 



It is not hard to see that VK/) an d are adjoints of each other (see for instance [2] for details). 

Furthermore, they satisfy the canonical (anti) commutation relations 

$(f),r(g)} T = (/, 5)1, $*{M*{g% = o, (A.i) 

where [A, B] T ■= AB =F BA, and ^ = or ip. In order to simplify notation, we usually identify cl 
with c, where c G C. 

For our purposes, it is more natural to work with the rescaled creation and annihilation operators 



:= —P 



1 

where N > 0. We also introduce the operator-valued distributions defined formally by 
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where 5 X is the delta function at x. The formal expression ip^{x) has a rigorous meaning as a densely 
defined sesquilinear form on T± (see [19] for details). In particular one has that 

= J dx f(x)Tp N (x) , 4>* N (f) = J dx f(x)$x(x) . 

Furthermore, the (anti)commutation relations (jA.ljl imply that 

$ N (x),r N (y)) T = jf8{x-y), $t(x)M(y% = 0, (A.2) 



B The Limit e — > in Loimria 1 6.7 



What remains is the justification of the equality in (|6.24p for e = 0. Our strategy is to show that both 
sides of (I6.25P with e > converge strongly to the same expression with e = 0. 

We first show the strong convergence of G[ k ' l ^' £ (a^). Let $ ( - n - ) £ and consider 

-Wy, || = H^i^i^-ijWije-^^WH < ||W-ye-^°$W||. 

Since the right side is in L 1 ([0,t]), we may use dominated convergence to conclude that 



Now 



lim fdsWiW^-Wij^^W = 0. 



ds I ds' \\Wt 3S W§ jls ,^ - Wij ta W w *V>\\ 
o Jo ' 

< f ds rds'ii^^^^w-^w-^^wii 

Jo Jo 

+ [ ds [ da / ||W^ ja Wi/ J v,,$W-Wi J> Wi/j/ ia /$(»)||. 
io io 

The first term is bounded by 

(^) 12 / ds 1l w w $(n) - w w* (n) ll - o, 

The integrand of the second term is bounded by 2|| Wy )S Wj/j' )S /3>( n ) || £ L 1 ([0, t] 2 ), so that dominated 
convergence implies that the second term vanishes in the limit e — ► 0. A straightforward generalization 
of this argument shows that 

G (M),^ a ( P ))$(p+fc-0 _^ G (M)( a (p)) $ ( P +fc-0 ) 
as claimed. Since the series (I6.24p converges uniformly in e, we find that 

oo k ^ oo k ^ 

EEf^(^(° w )) $(b) - EEfM^^))^. 

fc=0 Z=0 fc=0 Z=0 

as e — > 0. 
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Next, we show that e ltH N&( n ) — > e l4 ^jv<j)( n ). This follows from strong resolvent convergence of 
H £ N to .ff/v as e — > by Trotter's theorem [18]. Let W £ ■= Yli<j W £ j, and consider 

N \\(H £ N - i)- 1 ^ - (H N - i)- 1 ^^!! = \\{H £ N - i)-\W - W £ )(H N - i)" 1 ^!! 

< \\(w -w £ )(h n -iy^^W . 

Clearly < n ) := (H N - i)" 1 ^") is in the domain of H^. By the Kato-Rellich theorem [19] . is in 
the domain of Wij for all i,j. Therefore, 

\\{W l ,-W £ j ){H N -\)- l ^\\ = \\l {lWi3l>£ - 1} W l3 ^\\ 

as e — ► 0. Therefore 

e 1 ^ lv(a (p) ) e-" H ^$ (n) -» e 1 ^ A N (a {p) ) e~ itHN 
as e — > 0, and the proof is complete. 
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